A class of highly entangled many-body states that can be efficiently simulated 
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We describe a quantum circuit that produces a highly entangled state of N qubits from which one 
can efficiently compute expectation values of local observables. This construction yields a variational 
ansatz for quantum many-body states that can be regarded as a generalization of the multi-scale 
entanglement renormalization ansatz (MERA), and to which we refer as the branching MERA. In a 
lattice system in D dimensions, the scaling of entanglement of a region of size L D in the branching 
MERA is not subject to restrictions such as a boundary law L D_1 , but can be proportional to the 
size of the region, as we demonstrate numerically for D = 1, 2 dimensions. 
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Simulating quantum many-body systems with a clas- 
sical computer is generically hard, in that the computa- 
tional resources needed to describe the state of N quan- 
tum systems (e.g. N quantum spins) grow exponentially 
with N. However, in some appropriate limit, quantum 
systems effectively become the classical entities of our 
everyday world. There are several good reasons to inves- 
tigate the conditions under which quantum many-body 
systems can be efficiently simulated by classical means. 
At a fundamental level, the cost of classically simulating 
a many-body system may be used to define an opera- 
tional notion of quantumness or classicality. In addition, 
such studies shed light on where the power of quantum 
computation [1] lies, while they may guide the develop- 
ment of new quantum algorithms. At a more practical 
level, identifying classes of many-body states that can be 
efficiently simulated classically has led to a host of new 
tools [2|-[5| to investigate quantum many-body systems, 
which are of interest in a wide range of research areas in- 
cluding condensed matter physics, high energy physics, 
quantum chemistry, nanotechnology, etc. 

Examples of quantum evolutions that can be efficiently 
simulated include those that entangle only a constant 
number n of the N systems at a time [6[; and those 
that generate a restricted amount of entanglement be- 
tween any left-right bipartition of the N systems when 
ordered in one dimension [7[. Another important set of 
simulable quantum evolutions is obtained by generating 
it with a restricted, non-universal set of transformations 
(e.g. Clifford group [8j, and free boson and fermion evo- 
lutions, generated by quadratic operators [9[). Finally, a 
quantum circuit on TV two- level systems or qubits, can 
also be simulated classically, in the sense that expecta- 
tion values of local operators can be efficiently computed, 
by restricting the causal structure of the circuit, as is 
the case of the multi-scale entanglement renormalization 
ansatz (MERA) [Toj| , which is used as a variational ansatz 
for the ground states of many-body Hamiltonians. In the 
quantum circuit for the MERA, constraints on the causal 
structure result in a restricted amount of entanglement in 
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FIG. 1. (Color online) Quantum circuit that transforms 
the product state 10)®^ into an entangled state of N 
qubits. It contains 0(N \og(N)) generic gates (and 0(N 2 ) 
swap gates) organized in about T = log 2 (iV) horizontal lay- 
ers, and O(N) vertical branches. The circuit is self-similar 
and has a peculiar causal structure, such that the expecta- 
tion value of a local operator O can be efficiently 
computed. 



the many-body state. For instance, in a lattice in D > 2 
dimensions, the entanglement entropy S l d of a region 
of L D sites scales as and thus obeys a boundary 

law. However, the ground states of certain many-body 
systems display violations of the area law, such as the 
scaling L D ~ X logD 11, 12j], an d can therefore not be rep- 
resented by the MERA. 



Here we introduce a quantum circuit where, as in the 
MERA, the causal structure is restricted in such a way 
that one can efficiently compute expectation values for 
local operators. However, in contrast with the MERA, 
the scaling of entanglement entropy can reproduce viola- 
tions of the boundary law, all the way up to a bulk law, 
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FIG. 2. (Color online) (a) Computation of state |^t+i) from 
two copies of state \^ t ) (for t = 2), by means of a unitary 
transformation V made of swaps and two-qubit gates, (b) 
Computation of state |^t+i) from state |^t) (for t = 2), by 
means of an isometric transformation W . 



L D , which we demonstrate numerically for D = 1,2 di- 
mensions. We start by describing the circuit and proving 
its simulability in D = 1 dimensions, with the D > 1 
extension being straightforward. 

Quantum Circuit. — We consider a family of quantum 
circuits on N = 2 T (for a positive integer T) qubits or 
quantum wires, and 0(N\og(N)) generic unitary gates 
acting on two nearest-neighbor wires, supplemented with 
0(N 2 ) swap gates that exchange the position of two 
wires. Figure [T] shows the quantum circuit for T = 4. A 
discrete time t, t £ {0, 1, 2, • • • , T = \og 2 (N)}, organizes 
the generic two-qubit gates into (double) layers. The cir- 
cuit transforms the initial unentangled state 10)®^ of N 
qubits at time t = into an 7V-qubit entangled state 
at time t = T, 
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The quantum circuit has a branching, self-similar 
structure such that, for instance, the circuit for N qubits 
is obtained by merging two circuits for N/2 qubits, and 
so on. For notational simplicity, first we will describe the 
action of the circuit assuming that the gates on any two 
of its branches are the same (in Lemma 2 below, we will 
lift this assumption). At t = 0, the N wires are grouped 
into N/2 nearest neighbor pairs and, for each pair, the 
product state |^ ) ® |^o) (where \Sfr ) = |0)) is trans- 
formed into an entangled state of the two wires, so 
that applying the first row of gates has the net effect of 
transforming N copies of |^o) into N/2 copies of 
l^o)®^ —> \^i)® N ^ 2 . In general, at time t, a number 
of swaps and a (double) layer of two-qubit gates act on 
pairs of sets of 2 l wires, where each set of wires is in 
an entangled state |^t), to produce an entangled state 
of 2 t+1 wires, 

|* t ) |* t ) |* t+ i) = V(\* t ) ® |*t» , (2) 

see Fig. [2fa), so that from time t to time t + 1 the state 
of the N wires evolves as 




FIG. 3. (Color online) (a) The past causal cone of three out- 
going wires, indicated by the shaded region of the circuit, has 
the shape of a binary tree with 2 s branches at time t — T — s. 
In each of these branches, at any given time t it only involves 
three contiguous wires, (b) Accordingly, the reduced density 
matrix pr = p (here, T = 4), can be obtained by comput- 
ing a sequence of T — 1 reduced density matrices on three 
wires (here, {/>2, P3, Pa}), see Eq. [5] (c)-(d) p2 in terms of a 



small tensor network, (e) Computation of pt+i from pf 
applying gates and tracing out three wires. 
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Thus, the circuit produces the sequence of states 

|* 0> ®" _> W |* T _ l} ®3 _> |* T ), (4) 

where \^t) = |^) is the state of the N outgoing wires of 
the quantum circuit. For each value of t, state \^ t ) is a 
highly entangled state of 2 l qubits, see Fig. HJ 

Efficient computation of local expectation values. — A 
key property of the quantum circuit above is that the ex- 
pectation value (\I/|0|\I/) of any local operator O (which 
for convenience we assume to be supported on three con- 
tiguous outgoing wires of the circuit (l3[) can be com- 
puted efficiently. This is in sharp contrast with the 
same computation for a generic state of N qubits, which 
would have a cost exponential in N. This remarkable re- 
sult follows from the two lemmas below and the equality 
(\£|0|\£) = tr(pO), where p is the reduced density ma- 
trix for three contiguous outgoing wires of the quantum 
circuit with N(= 2 T ) wires, obtained by tracing out the 
other N — 3 outgoing wires. To discuss computational 
costs with more generality, we consider that each wire 
corresponds to a x-di mens i° na l quantum system, with 
X = 2 for qubits. 

Lemma 1 (particular case). — When any two branches 
of the circuit have the same gates, the reduced density 
matrix p can be obtained with computational time and 
space scaling as 0(x 12 log(N)) and 0(x 8 ), respectively. 

Proof. — Let us consider the past causal cone of three 
contiguous out-going wires, defined as the subset of wires 
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and gates that can influence the reduced density matrix 
p, see Fig [3(a). If we cut this causal cone into time 
slices, then at t = T we find the three out-going wires; 
at time t = T — 1 there are six wires; in general, at time 
t = T — s there are (at most) 3 x 2 s wires. We will 
obtain p by computing the state of these sets of wires, 
progressing from t = all the way to t = T. Importantly, 
the factorization of Eq. [4] for the state of all the wires 
implies that the density matrices inside the causal cone 
(for t > 2) factorize analogously as 

where p t is a reduced density matrix of three contiguous 
wires at time t, and pT = p, see Fig. [3] (b). Therefore 
p can be simply obtained by computing the sequence of 
T — 1 reduced density matrices {p2,p3, — m , Pt}- The 
first density matrix p2 is obtained from the entangled 
state |^2) of four wires by tracing out one of them, as ex- 
plained in Figs. 13(c)- (d). This can be seen to require time 
0(x 6 ) and memory 0(x 4 ). Then, for t = 2, • • • , T - 1, 
the density matrix is obtained from two copies of 
pt, as indicated in Fig. [3(e). This requires computa- 
tional time and memory that scale as 0(x 12 ) an d 0(x 8 ), 
respectively. Thus, the overall computational time and 
memory grow with x and N as 0(x 12 log(iV)) and 0(x 8 ) 
respectively, completing the proof. 

Lemma 2 (general case). — When all the gates in the 
circuit are allowed to be different, the reduced density 
matrix p can be obtained with computational time and 
space scaling as 0(x 12 N) and O(x 8 log(A0) Q- 

Proof.- The proof is analogous to that of Lemma 1. 
At time t the state of the N wires still factorizes as N/2 t 
states of 2* wires as in Eq. [4j but each set of 2* wires 
is now in a different state. Consequently, at time t the 
reduced density matrix inside the past causal cone still 
factorizes as N/2 f density matrices of three wires as in 
Eq. [5j but now each set of three wires is in a different 
mixed state. One such density matrix at time t + 1 is ob- 
tained by combining a pair of density matrices at time t, 
as in Fig. [3(b) (where branches are labeled by a sequence 
of A's and B's), and this has to be done for N/2 t+1 pairs. 
Adding the computational time for all values of t, we see 
that p can be computed with time and extra memory 
scaling as 0(x 12 N) and 0(x 8 log(iV)). 

Variational many-body ansatz. — Consider now a one- 
dimensional lattice C made of N sites, where each site 
corresponds to a qubit or, more generally, a %-level sys- 
tem, and a local Hamiltonian H = Hi that decom- 
poses as a sum of local terms Hi acting on at most a 
block of three contiguous sites including site I. A nat- 
ural application of the quantum circuit of Fig. [T] is to 
use the resulting state \^f) to define a variational class 
of states for the A^-site lattice £, which we could use, 
e.g. to approximate the ground state of Hamiltonian 
H. We refer to this ansatz as the branching MERA, 
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block length, L 

FIG. 4. (Color online) (a) Quantum circuit corresponding to 
the MERA. (b) Quantum circuit corresponding to a branching 
MERA with less branches than that it Fig. [lj (c)-(e) Holo- 
graphic trees for the circuits in Figs. 2(a)- (b) and Fig. [TJ The 
plot shows the scaling of entanglement entropy in a block of L 
sites evaluated numerically from a branching MERA created 
with random, free fermionic gates. Solid lines display a best 
fit over the indicated functional form. 

for reasons that will become clear shortly. The varia- 
tional parameters are contained in the gates of the quan- 
tum circuit (excluding the swaps), with each gate de- 
pending on 0(x 4 ) continuous parameters. In the general 
case of lemma 2 there there are 0(N log(iV)) independent 
gates [down to O(N) in the particular case of lemma 1]. 
Therefore, depends on 0(x 4 N\og(N)) [respectively, 
0(x 4z N)] variational parameters [15|- For any choice of 
these parameters, we can efficiently compute the energy, 
= (^f\H\^f) = Y,i (^\Hi\^), since we can efficiently 
compute each of the O(N) terms (^\Hi\^/). The varia- 
tional parameters in the branching MERA can then be 
optimized, e.g. by adapting the optimization scheme of 
Ref. [16|. 

The regular MERA, see Fig. HI (a), can be regarded as 

a particular case of the branching MERA, in which the 

gates for all but one branch are set to be the identity gate. 

Notice that in the regular MERA, at time t an entangled 

state \S&t) of 2 £ wires and the product state \0)® 2 of 2 l 

fresh wires are mapped into an entangled state of 
2 t+i 

wires, 

|*t> -> |*t+i> = W\9 t ), (6) 
where the isometric transformation W is obtained from 
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FIG. 5. (Color online) (a) Local detail of the branching 
MERA in D — 2 dimensions. Four different sub circuits are 
combined into one by means of a double layer of four-wire 
gates, (b)-(e) As in D = 1 dimensions, one may remove 
some of the branches by setting some of the incoming wires 
to a fixed state |0). This results in holographic trees with a 
branching factor ranging from one (regular MERA, (b)) to 
four (maximal branching, (e)). The plot shows the scaling 
of entanglement entropy in a block of L 2 sites evaluated nu- 
merically from a branching MERA created with random, free 
fermionic gates. Solid lines display a best fit over the indi- 
cated functional form. 



the unitary transformation V by fixing the state of 2 £ 
incoming wires to |0), W = VlO) 02 , see Fig. [2jb). In 
this case, a local expectation value (^\0\^) is also ob- 
tained by computing a sequence of reduced density ma- 
trices contained in the past causal cone of the support of 
the local operator O, but now the causal cone contains 
only one branch, and therefore Eq. [5] is replaced with the 
simplified sequence 



P2 -> P3 



pT-l 



(7) 



From a condensed matter perspective, one step of the 
regular MERA (from time t+1 to time t) can be regarded 
as implementing a coarse-graining (isometric) transfor- 
mation that maps the state |^t+i) of a lattice of 
2 £+1 sites into the state \^ t ) of a coarse-grained lattice 
of 2 t sites, and thus implements a real-space renormaliza- 
tion group transformation [5[. In contrast, one (reversed) 
step of the branching MERA can be regarded as imple- 
menting a decoupling (unitary) transformation that 
maps the state |^t+i) of a lattice of 2 £+1 sites into the 
decoupled state \^ t ) ® \^t) of two coarse-grained lattices 



of 2 £ sites each, thus disentangling two sets degrees of 
freedom that become decoupled at a given length/energy 
scale (l7j . 

More generally, setting a subset of branches to the 
identity, as in the example of Fig. H^b), produces other 
interesting subclasses of many-body states, all of which 
can be efficiently manipulated. These different subclasses 
of branching MERA can be roughly characterized by 
their holographic tree [17j , which represents their branch- 
ing structure, see Fig. H^c)-(e). For a generic choice of 
variational parameters, the holographic tree largely char- 
acterizes the structure of entanglement in the branching 
MERA. For instance, focusing on the entanglement en- 
tropy S l of a block of L contiguous sites (L < N/2), 
we find that different holographic trees imply different 
scaling of 5l, see Fig. |4j This scaling ranges from a 
boundary law, where Sl saturates to a constant, all the 
way up to a bulk law of the branching MERA of Fig. 
[TJ where Sl grows as « L; and it includes, as particular 
examples, the logarithmic scaling Sl ~ log(L) character- 
istic of the regular MERA in Fig. HJ^a), as well as a less 
orthodox scaling Sl ~ log 2 (L) of the branching MERA 
in Fig. \Mh). 

The generalization of the branching MERA to D > 2 
space dimensions is straightforward, see Fig. [5j and 
results in an ansatz where the scaling of the entangle- 
ment entropy of a block of L D sites also ranges, de- 
pending on the holographic tree, anywhere between a 
boundary law S l d ~ L D ~ X and a bulk law S l d « L D , 
including logarithmic violations of the boundary law, 
S l d ~ L D ~ X logL. The latter corresponds to the scaling 
of entanglement entropy in certain compressible gapless 
phases in D > 2 dimensions, including Fermi liquids [11 ] 
and spin-Bose metals 12], suggesting that the branch- 
ing MERA may offer a suitable representation for those 
phases. 
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